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1. INTRODUCTION 

Global warming and lowering fossil reserves have pushed many countries to move towards the use of 
renewable energies, which are considered clean energy [1]-[4]. In this context, the major areas of research in 
the world have been focused on the problem of optimizing wind energy captured by the turbine. Admittedly, 
the control of a wind turbine equipped with a fixed-speed alternator is very simple but can no longer ensure 
optimization of the aerodynamic energy captured. [5]-[8]. Because of this, variable speed operation was 
introduced. In this context, the MPPT control is called to maximize the wind energy extraction. Various MPPT 
control models were based on the standard model [9], [10] which neglects the nonlinearity of the magnetic 
characteristic as well as the iron losses [11]-[14]. 

For these reasons, this article develops the principles of the high performance MPPT control for DFIG 
based on the optimization of the currents passed through the stator and rotor. Optimality is grasped in the 
context of maximizing torque and decreasing circulating currents [15]-[19]. The theoretical results are justified 
by the simulation which confirms that the developed controller has a certain robustness to disturbances. For 
some authors, the reference rotor flux is given the nominal value (or by defluxing). In such a case, the 
electromagnetic torque is close to its nominal value and the efficiency of the machine is maximum. But, in 
reality, the electromagnetic torque is not always fixed in a precise value and can have strong variations 
(depending on the mechanical torque applied) [20]-[23]. So, if the aerodynamic torque is small, there is 
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untapped stored energy in the stator and rotor inductances reducing the energy efficiency of the machine [24]- 
[26]. The article is structured as follows: the model of the wind turbine taking into account the nonlinear 
magnetic characteristic as well as the iron losses are described in section 2. The formulation of the problem, the 
design of the reference signals, and the development of an improved MPPT controller based on backstepping 
technique which exploited the magnetic characteristic nonlinearity are established in section 3. The robustness of 
the control law is presented by simulation in section 4, a conclusion is illustrated in section 5. At the end of the 
article, the list of references and the bibliography are presented. 


2. WIND TURBINE MODEL IN PRESENCE OF NONLINEAR MAGNETIC CHARACTERISTIC 
AND IRON LOOSES 
2.1. Turbine model 
The aerodynamic wind power available at the blades is given by [19]: 


Prero= - prR? Cp (A, B)vw ( 1) 


with Bisa blade pitch angle which will remain constant for our case study, Cp(A) is the coefficient of 
performance (Cp(A, B) “& Cp(A)),R is the radius of the wind turbine rotor, p is the air density, vw is the wind 
velocity and A is the tip speed ratio, its expression is: 


(2) 


2.2. Nonlinear dynamic DFIG model 
2.2.1. Iron losses 

Generally, the iron losses of DFIG are most frequently modeled by of the equivalent iron loss 
resistance Rm associated in parallel with the magnetizing inductance X. although, in this paper, the proposed 
equivalent circuit of DFIG is inspired from [19], where the iron losses of DFIG are modeled by an equivalent 
resistance R,, associated at the same time in parallel with the stator cyclic leakage inductance L,, and mutual 
cyclic inductance X. In fact, the simplicity and precision of the calculations of this technique justifies this choice 
[13], because, if we use the theorem of thevenin, we will obtain an equivalent circuit like that of the standard 
model (without iron losses). Figures 1(a) and (b) show the proposed equivalent circuit of the DFIG in the 
reference frame aand B respectively using concordia transformation. So, the iron resistance is a variable 
parameter which was identified experimentally on a 4 kW DFIG driven by a DC motor. The procedure 
described in [19] was used to determine the values of the iron loss resistance as well as the magnetizing 
reactance [22]. The iron resistance characteristic is shown in Figure 2. 

From the Figure 1, the parameters Rs, Rm and Vsag will be replaced by an equivalent voltage source 
Vsrag in series with an equivalent resistance Rsr, from this, we are calculated in the stationary a- B reference 
frame the thevenin’s equivalents parameters as [23]: 


RsR, R R 
Rep = ; ViTa = ——— Vsa | VsTR = ——Vsp; 
sT Rg+Rm” sTa Rs+Rm sa > YsTB Rs+Rm sß > 
F Rs+Rm - Vsa.: Rs+Rm - Vsp 
ista = ia t; irg = isg + — 3 
sTa Rm 5% Rm” sTB Rm sp Rm ( ) 


(b) 


Figure 1. Proposed DFIG model in the (aß) reference frame: (a) a —axis and (b) B —axis (iron loss 
resistance associated in parallel with the stator cyclic leakage inductance and mutual cyclic inductance) 
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Figure 2. Experimental measurement of iron loss resistance equivalent characteristic 


2.2.2. Nonlinear magnetic characteristic 
We considered that the grid is stable and the stator flux aligned with the d axis. This gives: 


Vs 
Vsq = Vso Psq = 0 , and Ọsa = Ps * = (4) 


Ws 


Therefore, the stator and rotor flux can be rewritten as: 


Psa = Ps = Lsisra + Xtra (5) 
Psq = 0 = Lsisra + Xirqa (6) 
Pra = Lrira + Xisra = LrOira + (7) 
Prq = Lrirg + Xsrq = LrOirq (8) 


On the other hand, the stator and rotor currents contribute to the generation of the magnetizing flux along with 
the magnetizing current, denoted by: 


iu? ina = İsra + İra» ing = isra + irq (9) 
and we have Ppa = Xina: Pug = Xing (10) 


Where „is the magnetizing air-gap fluxes. In that case, the flux in (5), (6), (7) and (8) are given in the dq 
reference as [13]: 


Ps = lsolsra + Pua (11) 
Psq = 9 = Isoisrq + Pug (12) 
Pra = letra + Pua (13) 
Prq = lrotrg + Pug (14) 


The leakage inductances are grouped in the stator side, we pose lseq is the equivalent inductance to the stator 
side of both the stator and the rotor leakage brought, that gives l,o = 0 and lso = lseq. Next, the inductances 
L, and L, can be rewritten like that: 


Ls = lseq + X and lL, = lpo +X =X (15) 


Figure 3(a) give the variations of mutual inductance as a function of the corresponding magnetizing 
current. We suppose that the magnetizing flux direction remains unchanged and that it saturates only in 
amplitude, which allowed us to admit that ọ„and ipare in phase. Figure 3(b) shows the magnetic characteristic 
®, =A(1,,) where ®, and I, are the norms values. Of the above, the instantaneous values pand ipare given 
by: 
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Pu= T Íy (16) 


X= (17) 


First, we note some experimental points of the DFIG magnetic characteristic to make the approximation of the 
mutual inductance X(1,,). Then, we can easily deduce the magnetizing flux norm ®, (with ®, = ®,) in function 
of the magnetizing current I. 


Rotor flux norm 4, (Wb) 


0.04 i i i i i i 
1 2 3 4 5 6 7 8 re 
Magnetizing current Ip(A) Magnetizing current Ip(A) 
(a) (b) 


Figure 3. Experimental measurement of DFIG: (a) magnetizing inductance and 
(b) nonlinear magnetic characteristic 


2.2.3. Improved model of DFIG 
From the above and taking into account the iron losses and the nonlinearity of the magnetic saturation, 
the improved model of DFIG can be rewritten in the dq reference frame as follows: 


; i T 
f = f@wy) = [f,@wy) fœ uy) f œu y)] (18) 
h(x) = [hy (x) hE) hz (x)]" 
with 
X= [x1 X2 X3 X4 xs] = [isra isTq Pra Prq a” 
u= [Vsra VsTq ie y= [vra Vrq 7 
A . 1 1 1 
—Nisra + Wsisrqg + T'Prq + La, Pra + Lo, “STA = leg ra 
fi . ; 1 í 4 X 
h —Uslsra — NlsTq — beg Pra +I Prq + Deg STA = leg "a 
‘5 = Rrisra = lseq! Pra T IVsPrq + Urq (19) 
4 RrisTq = GDsPra — lseq! Prq F Vrq 
fs H . . f Tmec 
T (Psaira — Psalrq )-72 E 
2 
Where: lseg =oL, ; 0 =1- 3 ppe eT pe, u= sea . With: J : total inertia constant ; 
LrLs lseq lseq% lseq T 


f: viscous friction coefficient; (Ps, g,): stator , rotor flux components; 2: DFIG rotor speed;(R, , R»): stator 
and rotor resistances. 
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Figure 4. Proposed DFIG equivalent circuit accounting the iron losses and the nonlinearity of the 
magnetic saturation: (a) d axis and (b) q axis 


Figures 4(a) and 4(b) show the proposed circuit of the DFIG taking into account the iron losses and 
the nonlinearity of the magnetic saturation in the reference frame dandq respectively using Parck 
transformation. In the above model, the parameter T is the unique parameter that varies in function of magnetic 
state of the DFIG, this relation has been given a polynomial approximation [24]-[26], i.e.: 


T=K(®,) , with : K(®,)=ao+ ai®, +a,” +....4 an,” 


The ais are constant coefficients, identified with polyfit MATLAB function. So, the magnetic 
characteristic of the DFIG is denoted ®, = y(1,,) where I, denotes the norm of the magnetizing current and 
P, = P. From the experimental points (1,,,®,), we deduce a polynomial function approximation of the 
characteristic (®,, I) as depicted Figure 5. The parameter values of DFIG used in experimental tests are given 
in Table 1. 
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Figure 5. The characteristic (®,, T) 


Table 1. Parameters of DFIG 4 kW 


Characteristic Value 
Stator resistance 1.250 
Rotor resistance 0.172 
Mutual inductance 0.0772 H 


Rotor leakage inductance 0.0018 H 
Stator leakage inductance 0.00096 H 
Generator inertia 0.33 kg m? 
Viscous friction coefficient _0.001 m/s 


3. REFERENCE ROTOR FLUX CONSTRUCTION AND THE ENHANCED OPTIMAL 
CONTROL 

The present section presents a new and optimal approach that properly accounts for the nonlinear 

magnetic characteristic to control DFIG. Notice that in previous works, the standard model can be referred to 

linear model (LM), from the nonlinear model (NLM), and this while keeping the parameter I constant. Thus, 

we describe a new control strategy based on a method for optimizing the rotor flux reference to minimize stator 
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and rotor current absorbed [23]-[26]. the performance of this optimization law is validated by means of 
simulation. 


3.1. Problem formulation 

The problem of controlling variable speed wind turbines is addressed. The majority of previous works 
considered that the reference of the reactive power of the stator was null and the magnetic characteristic was 
linear. However, the machine operation mode is not optimal of the large electromagnetic torque with Qs=0.In 
the present section, we will develop a new optimal MPPT control strategy taking into account the nonlinearities 
of the machine (in order to lower the total currents consumption in both the stator and rotor). 

The control objective is to follow the maximum of the wind power captured and to minimize the 
currents transited by generating an optimal rotor flux in line. In this context, the rotor flux norm ®, must track 
perfectly the flux reference Ọ -ref = H (ly Is) , where (1,,1,) are the rotor and stator currents norms and the 
function H will be resolute in order that P, = P -ref involve the minimization of the stator and rotor currents 
circulated in the induction machine for getting a maximum electromagnetic torque [23]-[26]. 


3.2. Optimal rotor flux reference 
Firstly, we find a relation between the rotor flux, stator and rotor currents norms. In this respect, one 
considers the DFIG model developed in (19); thus, the stator and rotor currents can be rewritten by: 


é L = : Ko 

id = sPrd-XPs — Prd _ Vs gS rq (20) 
oLsLy oX lseq Xo 

; — Ps- Yrd ; = X. 

lta = leq? lstq = Ta (21) 


On the other hand, the rotor flux can be rewritten as: 
Přa=X? (igra + ifa + 2israira) » rg=X° (iSra + leg + 2istgirg) (22) 


Then, using (20) and (21), we'll have: 


2 2 2 2 2 

. . PsPrd Pra Ps PsPrd 1 1 Pr Prq Psd 

istairg = ot — fra SE p SS = (a+ - +— 1 fsa (23 
shore OLylseq OLylseq eq eg PsPra eg OLylseq OLylseq OLylseq eq ( ) 

2 

i š Prq 

istgirg = — 24 
sTq*rq oLrlseq ( ) 

Accordingly, one has: 
2 es 2 2 2 . . . . 
p; = X (I; F I; F 2 (istara + istgirg) (25) 


We replace isaira , and isgirg by their expressions according to (23)-(24), we are the expressions of rotor flux 
norm in function to the stator and rotor currents norms: 


1 1 2 5 
of = UE + 17) +287 (pH SF) sara + sara) — 28° GS 2 oo 


oLylseq OLylseq 
The electromagnetic torque Tem can be written by: 


_ P(@saPrq-PsqPrd) 
Le (27) 


Then, (26) can be simplified to: 


2 
s 


1 2 2 IseqVs 
($+ )@- +) +2-@- es 


x? OXlseq Beg 


Z ) =0 (28) 


OLrlseq 
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Figure 6 shows the curves representing the sum of the squares of the current modules (12 + I?) versus 
the flux ®, , for various values of the electromagnetic torque Tem. At present, it is easy to observe that, for 
each value of the sum (12 + IŻ), there are several different operating points by the value of rotor flux norm ®, 
and Tem. For example, a torque Tem = 5Nm can be produced, for the generator 4kw (Table 1), with: 

— A current I + I? = 490 A and a flux ©, = 0.2 Wb; 
— A current I2 + I? = 20 A and a flux ©, = 0.40 Wb; 

The best operation point from an energetic viewpoint, is the one which corresponds to the minimum 
of currents circulations. From Figure 6, it is usually seen that there is a unique couple (®,,(I2 + I?)) for any 
Tem, Which engages the least possible values of the stator and rotor currents. For our case study ® -ref = 
0.4Wb. 

The stator active and reactive powers are given by: 


Ps = Vora İsra + Vsrq İsTq (29) 
and 

Qs = Verqisra — Vera İsTq (30) 
For Vsrq = Vs and Vra = 0, (29) and (30) become: 

B = Vista GD 
and 

Qs = Visra (32) 


Introducing (20) and (21) into (31) and (32), we obtain: 


_ VX . 
Py = -T inq (33) 
and 
Qs = Ņ -Vy ia (34) 
Ls Ls 


Indeed using (33)-(34) and (20)-(21), we can deduce the expression of the rotor flux magnitude reference 


(Dre = | Pra + P7) according to Qs ref (stator reactive power reference) 


rér = [vs = 00, sep] + [seer] (35) 


The stator flux @, is assumed constant (the grid is considered stable). 


avec 


Is? + Ir? 


fi i i n i mone 
0.1 0.15 0.2 0.25 0.3 0.35 0.4 0.45 0.5 
Rotor flux norm o, (Wb) 


Figure 6. Sum of the squares of the current modules (12 + I?) as versus rotor flux ®, for different Tem values 
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3.3. Improved MPPT control 

The main objective of the control strategy is to develop an optimal control law in order to optimize 
the different operations points of the wind turbine [13]. For this, an enhanced backstepping controller has been 
designed to force the system to track their references ( Qref , Pr-ref) Where Nhref is the optimal rotor speed 
reference illustrated in Figure 7 and ®,_,¢¢ is the rotor flux reference (35). 


Ores (rad/s) 
50 


:— MPPT 


Windspeed 
0 2 4 6 8 10 12 (m/s) 


Figure 7. Optimal rotor speed reference in function the speed wind [13] 
For reasons of simplification (backstepping control of the rotor flux in only one step), the control of 


the rotor flux will be assured by controlling the stator reactive power described by (34). Firstly, we consider 
the following errors: 


e1 = Qef — 0 (36) 
e3 = Qs ref =Q; (37) 


Step 1. It follows from (19), (33), (34), (36) and (37) that the dynamics errors e, and e}are: 


, ; ; : u , f T, 
êi = ref — Q = Aref +5 Psairq Haee (38) 
? . ‘ . Rs d., Xd. 

e3 = Qs ref - Q; = Qs ref +V Tews ae 874 + Vs 7 aeird (39) 
r 5 XR; d., D 

e3 = Qs ref z s Za, mra Vs 7 ap ira (40) 


Indeed, considering the function Lyapunov candidate: 


1 1 
Vi =5et +363 (41) 
We deduce that: 
V,=e,é, + e3@, = —k et — ke2 (42) 


Where k, and k; are positive real parameters to be chosen as described in [13]. On the other hand, in (38), 
Psaltrq is considered as a virtual control signal which will temporarily be like the real control signal, the error 
in (38) assures the global asymptotic stability, letting: 


g aef J 5 f Tmec 
(Psaira) et = With: y, = al Ores = 72? + 7 = ke1) (43) 


As the quantity Psqirq can’t be equal to p, because isn't actual control signal. Indeed, we take the stabilizing 
function p; in the first. Step 2. Now introduce the new error: 


€2 = Hi — Psa irq (44) 
Then, the dynamic of the error e,are given as: 


é = —kye, — Tez (45) 
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So, the time derivative of the function V,is: 
V, = -kef — kze — Tee (46) 


The actual control signals V,q and Vg are choosing so that all errors (e4, €2, e3 ) converge to zero: 


: . d., 
€2 = W1 — Psa tra (47) 
: J A fA Tmec . d., 

and é, = F. (—Qref — 7? ET k1€1) — Psa z; İra (48) 
From (3) and (4), we can write: 

dad, 1 P X ds P 

uri = Lro (Vra _ Rrira Tx a T L,OWrirq) (49) 

ae 1 r j X 

utra = Lro Veg = R, rg — LO, ipa — 1, Or Ps) (50) 


Introducing (25) and (26) in (27) and (24), we obtain: 


; J[_ ä f anih o ; F X 
e2 = u (Gree ~ 7? F F F kže, ) F ke2 _ ee Vig E Rrirq 7 L,OWrlyrq = 1, Ors) 


This can be condensed as: 


: V; 
bap (51) 


With u = Sires = Là + T + kfe,) + k,e, ~ esd (—Rrirq — LOW; tpg — * Ps) 
T S 


ee _ XRsVrq XRsVs ( . . x ) XVs 
and e3 = Qs ref Vy Ls?Lrows | Le2Lyows R,irg + L, 0, irq + Ls W, Qs) + LsLro Vra + 
XVs . X d@s Fi 
Bie. Rripa Ls dt + LOW, irg ) (52) 


Then, let us consider the Lyapunov function candidate as: 
1 1 1 
Vz = Set +7e%+7e3 (53) 


The control signals V,, and V,q can be deduced from (51) and (52) satisfying V, <0: 


Lro 
Yq = P (k2e2+42) (54) 
LsLro a XRsVy . X . X dos s 
V,a = ra |-kses — Qs ref + Vigo — Owrira + +o 0s)| + Rrira + ide Ly OWrirg (55) 


Where k, is a new positive constant. In fact by substituting (54) and (55) in (53), the time derivative of (53): 


Which ensures the global asymptotic stability of the system. 


4. SIMULATIONS RESULTS 

In this section, we will present the importance of the new control strategy (non-linear model with 
optimal flux NLMOP), introducing the optimal rotor flux reference illustrated by Figure 8, compared to the 
control strategies using the constant rotor flux reference (case of the linear model with constant flux LMCF or 
nonlinear model with constant flux NLMCF). The comparison is performed using DFIG and Turbine 
parameters considered in Tables | and 2. 
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Rotor flux reference $ |, (Wb) 


Electromagnetic torque Tem (Nm) 


Figure 8. Characteristic of the rotor flux reference obtained by interpolating the points (Tem, Pr-ref) 


Table 2. Turbine parameters 


Characteristic Value 
Turbine radius, R 3m 
Gear box ratio 1:8.5 


Rated wind speed 10.5 m/s 
Number of blades 3 

Power 5S kW 
Rotor inertia 0.05 kg m? 
Cut-in wind speed 4.25 m/s 


4.1. Hegemony of the NLMOF strategy over the LMCF strategy 


The simulation protocol is designed in order to test the performance of the HP-MPPT controller (used 
the NLMOF) operating under different conditions by considering two steps of the wind speed (Figure 9) and load 
torque (Figure 10). Recall that the NLMOF regulator involves the optimal flux reference depicted in Figure 8. The 
HP-MPPT control laws defined by (54) and (55) receives the following design parameter values, which have 
proven to be practical: kı = 1000, kz = 300 and k = 50. 

The optimal specific speed 8,1 with MPPT control are presented in Figure 11. The dynamic responses 
of the optimal flux NLMOF and constant flux LMCF regulator are presented in Figure 12. The control of rotor 
flux with the NLMOF and LMCF strategies are presented in Figures 13 and 14 respectively. According to the 
results obtained in Figure 15, it is confirmed that under all operating conditions, the NLMOF regulator requires 
smaller current than the constant flux reference regulator LMCF. Figure 16 confirms the performance of our 
regulator by the pace of joule effect losses for LMCF and NLMOF strategies. The active stator and rotor power 
are presented in Figures 17 and 18 respectively for LMCF and NLMOF regulators. The dynamic and the 
performance under different operating parameter changes of the proposed regulator have been affirmed further 
by Figures 19 and 20 which presented the reactive stator and rotor power respectively for LMCF and NUMOF 
controllers. 


Wind speed (ms) 


mechanical torque (Nm 


Time (s) 


Time (s) 


Figure 9. Wind speed Figure 10. Load torque applied 


Int J Pow Elec & Dri Syst, Vol. 13, No. 2, June 2022: 1195-1208 


Int J Pow Elec & Dri Syst 


Time (s) 


Figure 11. Optimal specific speed with MPPT 
control 


==. Rotor flux reference |] 
Rotor flux mesured | | 


o ss Ae ea es a LS 
1 2 3 4 5 6 T 8 9 10 
Time (s) 


Figure 13. Regulation of rotor flux with the 
NLMOF strategy (The online flow reference 
generated (Tem, Pr-ref)). 
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Figure 15. Sum of squares of norm (12 + I?) for 
LMCF and NLMOF regulators 
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Figure 17. Active stator power (W) for LMCF and 


NLMOF regulators 
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Speed (rds) 
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Figure 12. The Speed reference (to ensure MPPT), 
the responses of the optimal flux NLMOF and 
constant flux LMCF regulator 


Rotor Flux (Wb) 


= Rotor flux reference 
Rotor flux mesured | | 


a an) 
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Figure 14. Regulating the rotor flux with the 
LMCFE control strategy (which imposes a constant 
flux reference) 


= NLMOF strategy 
=== = LMCFstrategy 


Joule loss (W) 


M00 E NE i 


Figure 16. Joule effect losses for LMCF and 
NLMOF regulators 
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Figure 18. Active rotor power (W) for LMCF and 


NLMOF regulators 
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NLMOF strategy 
-  LMCFstrategy 


—— NLMOFstrategy 
===- LMCFstrategy 


Stator Reactive power (VAR) 
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è 
S 
s 
Rotor Reactive power 
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Figure 19. Reactive stator power (VAR) for LMCF 
and NLMOF regulators 


Figure 20. Reactive rotor power (VAR) for LMCF 
and NLMOF regulators 


4.2. Hegemony of the NLMOF strategy over the LMOF strategy 

The simulation protocol is similar to part 4.1 (Figures 9 and 10). Recall that the NLMOF regulator is 
defined for (19) in combination with the optimal flux generator (shown in Figure 8). So to have the LMOF 
regulator, it is sufficient to make the following changes in (19): 
— The parameter T in (19) is held constant, equal toà r ~ 460 QH7?. 
— The rotor flux reference is given by the signal depicted in Figure 8. 
— The controller parameters are set to the following values: k; = 1000, k, = 300 and k; = 50. 
Figure 21 shows that the Sum of squares norm (12 + I?) is lower with the NLMOF regulator. Figure 22 shows 
the quality of flux tracking obtained with the NUMOF and LMOF control strategies. So, Figure 23 presents the 
active stator power for LMOF and NLMOF regulators. Also, from the energy point of view, is also illustrated 


by Figure 24 when the rotor power produced by NLMOF regulator is more than LMOF regulator. 
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Figure 21. Sum of squares norm (12 + I?) for 
LMOF and NLMOF regulators 
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Figure 22. Regulating of rotor flux for LMFO and 
NLMOF regulators 
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Figure 23. Active stator power (W) for LMOF and 
NLMOF regulators 


Figure 24. Active rotor power (W) for LMOF and 
NLMOF regulators 
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5. CONCLUSION 

The problem of controlling the DFIG in the presence of nonlinear saturation magnetic state are 
evoked. A new speed and rotor flux regulator has been developed using the optimal backstepping approach to 
optimize the currents circulated in the machine. Also, the proposed regulator comprises a reference rotor flux 
generator dependent on the nonlinear magnetic state. In addition, it is demonstrated by simulation that, under 
all operating conditions, that the current transited is really lower with the regulator NUMOF. The performance 
of the proposed control laws has been demonstrated by simulation studies performed in the 
MATLAB/Simulink® environment. 
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